Pair — correlations and the survival of superconductivity in and around 

a super-conducting impurity 
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The problem of the survival of superconductivity in a small super-conducting grain placed in a 
metal substrate is addressed. For this aim the pair correlations and super-conducting gap around 
and inside a negative-U impurity in one and two dimensions is calculated, in a discrete tight-binding 
model and a continuous model. Using a mean-field decomposition, it is shown that finite pairing in 
the grain develops when the system has a degeneracy between successive number of electron pairs, 
and thus may oscillate as a function of the chemical potential. For finite pairing in the island, 
pair correlations in the normal part exhibit a cross-over from being long-ranged to exponentially 
decaying, depending on the strength of interaction in the grain. It is shown analytically that there is 
a minimal island size under-which pairing vanishes which is different than that given by Anderson's 
criterion, and that it scales as a power-law with island size, rather then exponentially as in isolated 
grains. 

PACS numbers: 74.81. Bd, 73.63. Bd, 74.50.+r, 81.07.Bc 



I. INTRODUCTION 

While superconductivity (SC) on the nano-scale has 
been a long-standing issue, dating back to the seminal 
work of Anderson^, only in the last decade has technolog- 
ical advancement enabled the realization of such systems 
in experiment^. Since then, manifestations of SC on the 
nano-meter scale has been observed not only in isolated 
grains^ or granule on insulating substrate^, but also in 
inhomogeneous SC thin films, where well-separated SC 
and normal regions have been observed?. In such hy- 
brid systems, the proximity between the SC and normal 
phases gives rise to novel effects, mainly manifested in 
the local density of states (LDOS), which may be directly 
measured using scanning tunneling microscopy^. 

Encouraged by the technological advance, we ask the 
following question : what would be the properties of an 
ultra-small SC grain placed on a metallic (or a doped 
semi-conducting) substrate ? In such a case, would the 
gap in the grain still obey Anderson's criteria 1 , or will the 
proximity effect yield a new criteria for the destruction 
of SC in such a grain, as seen in, e.g. thin SC layers at- 
tached to normal layeri^ ? Furthermore, in such systems 
one expects that the SC properties of the grain would be 
strongly affected by the properties of the surrounding 
metal, and that the proximity to a SC grain would gen- 
erate pair correlations that would impinge on the local 
properties of the metallic area, such as its LDOS 9 . The 
above question is also interesting form a technological 
point of view, as Josephson arrays fabricated on metallic 
or semi-conducting substrate seem to have large techno- 
logical potential as nano-electrical devices. 

In order to examine these issues, a minimal model 
of a single SC grain placed in a clean metal matrix is 
studied, by means of a negative-U Hubbard Hamiltonian 
in which the interactions are confined to a small region 
in space (so-called "negative-U impurity"). Applying a 



Hartree-Fock-Gorkov mean-field decomposition^ leads 
to the Bogoliubov-De-Gennes (BdG) Hamiltonian 11 , 
which serves as a starting point in the calculation. 

The effect of the proximity between the SC grain 
and normal area is investigated by numerically solving 
a tight-binding BdG Hamiltonian in one and two dimen- 
sions. It is found that the pairing in the grain is strongly 
affected by the chemical potential (i.e. density) of the 
substrate, and that on the normal area pair-correlations 
may either be suppressed exponentially away from the 
impurity or be long-ranged, depending on the value of 
the attractive interaction in the grain. 

The dependence of the gap in the grain on its size is 
studied using a continuous version of the BdG Hamil- 
tonian. Solved analytically, the dependence of the gap 
on island size is found to diminish as a power-law rather 
than exponential (as in an isolated grair>i£), and the min- 
imal island size under-which SC vanishes in the graini^ 
is evaluated, and is found to depend on the properties of 
the substrate. 



II. THE NEGATIVE-U IMPURITY IN THE 
TIGHT-BINDING MODEL 

Let us start by examining the discrete tight binding 
model for the negative-U impurity. The model Hamilto- 
nian is 



(1) 



where t is the hopping clement, \x is the chemical po- 
tential and U > is the attractive interaction, which is 
only present in a single site at the origin (the negative- 
U impurity). By applying the Hartree-Fock-Gorkov 
decomposition 10 the BdG mean-field Hamiltonianii is 
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H.c. 



obtained, 

(2) 

where A = — [/(cL-cL ) is the pairing potential and 

&i = <5io2 ff I^KcLcot)/ 2 = <5io|C^|("-o)/2 is the Hartree 
shift. This mean-field approach is justified by noting that 
little is known about this system, and hence a prelimi- 
nary mean-field treatment is in place. The importance of 
maintaining the Hartree shift term, which naturally ap- 
pears from the derivation of the mean-field Hamiltonian, 
has been discussed and demonstrated in Reffl2l. where a 
similar decomposition was used to study the properties 
of a disordered superconducting sample. 

Introducing a Bogoliubov transformation one obtains 
from the Hamiltonian of Eq. @ the BdG equations^ for 
the quasi-particle (QP) and quasi-hole excitations 

v(ri), 
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(3) 

In Eq. © |u fe (r;) = -tJ2$ u k (ri+S) where 6 = ±x,±y 
and similarly for Vk(?i), and the energies are the QP ex- 
citation energies Ep. > 0. The pairing potential A and 
the electron density per site are to be determined self- 
consistently in terms of the QP amplitudes u(ri) and 
«(r 4 ) , 



A = |C/|5> fe (0K(0), 



2j> fc (r 4 )| 5 



(4) 



The pairing amplitude A is finite only on the negative- 
U impurity. However, the proximity to the impurity in- 
duces pair-correlations fa = {c\^c\^} even for i ^ 0, that 
is outside the impurity. 



III. RESULTS IN ONE AND TWO DIMENSIONS 

The pair-correlations are investigated by solving 
Eq. © numerically and self-consistently for a ID and 
2D metallic substrate, from which both the densities and 
pair correlations may be calculated. In Fig. [1] the pair- 
ing amplitude fo (bright points in Fig. []} is plotted as a 
function of [i and U for ID lattice of size L = 49 and a 
2D lattice of size 7x7 (Fig. HJa) and (b), respectively). 
In the calculation, hard-wall boundary conditions were 
taken, but similar calculation using periodic boundary 
conditions showed no qualitative change in the results. In 
ID, finite pairing is only visible along the line of constant 
density (n) « 0. In 2D stripes of finite pairing appear on 
lines of constant density, which are the boundary lines 
between sequential even (mean) number of electrons, JV, 
in the system. It is clear that the phase-space available 
for finite pairing is much larger in 2D than in ID. 
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FIG. 1: Color-plot of the pairing amplitude in the impurity, fo, as 
a function of chemical potential fi and interaction strength U for 
(a) one dimension and (b) two dimensions. Bright area indicates 
large pairing amplitude. 



In Fig. [5] the number of electrons N and the pairing 
potential A are plotted as a function of /i for a 13 x 13 
system with U/t = 3.5. As seen, the pairing potential 
(stars, right axes) is finite only at the transition between 
plateaus of constant even electron number N (diamonds, 
left axes). 
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FIG. 2: The pairing potential A (stars) and mean total electron 
number N as a function of chemical potential for a 13 X 13 system 
with U/t = 3.5. The pairing potential (stars, right axes) is finite 
only at the transition between plateues of constant even electron 
number N (diamonds, left axes) 



The reason for this behaviour of the pairing is that the 
Hamiltonian of Eq. ([2|) couples between states with no 
electrons and two electrons at the impurity. Thus, the 
self-consistent pairing is finite only when states with N 
and N + 2 electrons are degenerate at the Fermi energy. 
In ID this only happens when the occupation changes 
from N = to N = 2. In 2D, however, this degeneracy 
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is much more common, resulting in many regions of finite 
pairing in fi— U phase space, and hence in the oscillatory 
behaviour shown in Fig. [21 We note that while for 3D 
the computation is numerically demanding and will not 
be presented here, we expect similar behavior, with even 
more phase space available for SC in the grain. Such a 
case may be more relevant from the experimental side. 

Further insight may be gained by studying the depen- 
dence of A on the interaction strength U. Although from 
Eq. ([4J it would seem that the two are linearly dependent, 
this is not the case. Due to the Hartree term, U affects 
both the occupation of the impurity and the energy lev- 
els of the system, pushing the system in and out of the 
N — N + 2 degeneracy required for finite pairing. This is 
demonstrated in Fig. [31 where A (stars, right axes) and 
N (diamonds, left axes) are plotted as a function of U for 
a 9 x 9 system at fi/t = —2.4. This chemical potential 
corresponds to a low (« 0.25) electron filling, which is 
relevant for a semi-conducting substrate. A(U) is a non- 
linear (and non-monotonic) function, only finite above a 
certain critical interaction U c (n), in the transition region 
of N from N — 18 to N — 20. We note that by changing 
\x one may find finite pairing at higher electron densities. 
However, this is unlikely to occur above half-filling, as 
the Hartree term will suppress the pair function in that 
case. 
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FIG. 3: The pairing potential A(0) and mean total electron num- 
ber TV as a function of interaction strength U for a 9 X 9 lattice with 
fi/t = —2.4, showing that A(U) is a non-monotonic function, only 
finite above U c when TV fluctuates between TV = 18 and TV = 20. 
Inset: the numerical data (stars) in the region U ~ U c is fitted with 
the function A(U) <x (U/U c - l) x e(U - U c ) (solid line), yielding 
x ~ 1/2. 

In the region where the pairing in the SC island is 
finite, the proximity effect should yield pair correlations 
away from the impurity, /, 
fi on a chain of length L 
for two values of interaction, U/t — 1.5 (stars), within 
the energy band, and U/t — 2.5 (squares), outside the 
band. The chemical potential is adjusted for each value 
of interaction in order to maintain finite pair correlation 
in the impurity. When U lies outside the band, it is found 
that the hole excitations v(n) become localized, resulting 
in an exponential decay of the pair correlations. On the 



= (ct.^). In Fig. H we plot 
= 51— (normalized to unity) 



other hand, If U lies within the band the hole excitations 
are periodic, and generate long-range pair correlations. 
While this effect may be due to the finite size of the 
normal system, these long-ranged correlations may have 
a crucial effect on the global behaviour of a system with 
many negative-U impurities (so-called dilute negative-U 
model— i 17 ' 18 ), as they determine the effective Josephson 
coupling between the different impurities. 
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FIG. 4: Spatial structure of the pair correlations fi = (ct^ctj) 
in a linear chain of length L = 51 for two values of interaction, 
U/t = 1.5 (stars) and U/t = 2.5 squares. When U lies within 
the band there are long-range pair-correlations, but an exponential 
decay of the correlations for U outside the energy band. 



IV. TWO DIMENSIONS - CONTINUOUS 
MODEL 

The starting point for the following calculation is the 
continuous negative-U Hamiltonian (with H = 1) 

H = Y^ /d 2 r^(r)(-J-V 2 - M )^(r)+W/ , (5) 
where /i is the chemical potential, and 

Hi = j tPrCrCr^rJttJCr)*^)^) , (6) 

where U[r) < is a short-range attractive electron- 
electron interaction. The existence of a negative-U im- 
purity is modelled by limiting the interaction to a finite 
island-like region in spaced, i.e. 



U(r) 



\U\, rel 
0, else 



(7) 



where I is a disk with radius a around the origin. 

Using the BdG mean-field decomposition^, we substi- 
tute Hi by 



H A = I d 2 r(A(r)* t T (r)* t i (r) +H.c. 



(8) 
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where A(r) = |J/|(\&i-(r)\&|(r) is the pairing potential. 
The Hamiltonian TIa is now expanded with a being a 
small parameter (specifically kpa << 1), 



H A w {Aa 2 *}(0)*j(0) + H. 



(9) 



TyAW 2 (*}(r)*{(r) 



H.c. 



r=0 



where A = A(0). Since V 2 f/'j(r)-0{(r) |r=o ~ 
fe^^|(0)*J(0), The second term in Eq.l|IS]l is smaller 
than the first term by a factor of (fc/a) 2 and can be ne- 
glected, and we are left with the first term in Eq.(|10p 
which can be treated as a local perturbation. Notice 
that A is not assumed to be small, as it might not be. 

For typical metals, the condition /c^a << 1 yields an 
island of size of a few nano-meters, for which the contin- 
uum formulation is inappropriate. For semi-conductors, 
where the Fermi wave-length may be a few orders of mag- 
nitude larger than in metals, the continuum limit will 
still be valid. In semi-conductors which exhibit true su- 
perconductivity, due to the low carrier density the key 
role is played by inter-valley processes^. In the above 
model, on the other hand, attractive interactions take 
place only on the impurity and the sample as a whole 
need not become superconducting (on the contrary, it 
is assumed that it remains normal). Thus, local super- 
onducting correlations will probably not be detected by 
conventional transport measurement. However, the man- 
ifestation of pairing correlations can still be detected via 
local measurements of, e.g. the LDOS, in which a mini- 
gap should appear (see Eq. (fT2")) below). 

As a first step let us calculate the single particle 
LDOS in the presence of the impurity, given by p(r, 0) = 
— i3G r (r, r, 0) where G r (r, r', 0) is the retarded Green's 
function, which obeys the Dyson equation (depicted in 
Fig. El) 

G r (r,r',0) = g r (r,r',0) + A 2 aV(r,O,0) x 

x/(O,O,-0)G r (O,r',0) , (10) 
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FIG. 5: (a) Dyson's equation for the single particle Green's func- 
tion, (b) Dyson's equation for the anomalous Green's function. 



Where g r (r, r', 0) = -infl J (y / 2m(^i + 0) | r - r' |) is 
the Green's function for a clean 2D system, Jo is the zero 
order Bessel function and f2 is the 2D DOS at the Fermi 
energy. The solution for Eq. (fit)]) is 



G r (r,r',0) = ,f(r,r',0)- 



7 



Tltt 1 +7 



<f(r, O,0)<f(Q,r',0) , 
(11) 



where 7 = (irtlAa 2 ) 2 . This yields for the DOS 



p(r,0)~fi 1- 



1 + 7 



J o ( v / 2m(Ai + 0) I r |) ; 



(12) 



which shows a decrease (so called mini-gap) in the LDOS 
at the location of the island, and Friedel-like oscillations 
fro m the island. The oscillations persist on a length-scale 
~ y/n/li, independent of the pairing potential in the is- 
land, which only affects the depth of the minigap. Due 
to these oscillations in the LDOS, one expects that the 
charge density will re-distribute (and also exhibit oscilla- 
tions) around the SC-semi-conductor interface. One also 
expects that the density will re-distribute within the SC 
island. However, this effect cannot be probed within the 
point-island limit (Eq. (jTOJ) ) , and will be discussed in a 
future study. 

Next we turn to a self consistent calculation of A. 
Using the Dyson equation for the anomalous Green's 
function F(x, x', t) = (Tip\(x,0)i/}\(x',t)) (Fig.^b)) one 
finds 



F(x,x',0) 



Aa 2 
1 + 7 



g(x,0, 0)9(0, x', -0) . (13) 



The asymptotic properties of Bessel functions yields 
F(x) ~ 1/x far from the impurity, independent from the 
value of interaction strength, as expected in a N-SC in- 
terface. This is in contrast with the exponential decay in 
the tight-binding model. 

The self-consistency equation for the pairing potential 
reads (see, e.g., Ref. [ill) 

A(x) =-\U\ d0F(x,x,0) , (14) 

J —0D 

where 0b is the frequency cut-off of the interaction. Sub- 
stituting Eq. (fT3)) yields 



A(x) 



\U\Aa 2 



7 



1 + 7 



r)(x) 



(15) 



where rj(x) — J_^ D d0g(x, 0, 0)5(0, x, — 0). At the center 
of the impurity 77(0) = — 2(tt£1) 2 0d- Inserting this into 
Eq. (fT5|) results in an algebraic equation for A(0), 



1 = 2(7rft)V|[/|0 D 



(7rOa 2 ) 2 A(0) : 



1 + (7rf7a 2 ) 2 A(0) 2 
which is easily solved to give 

1 ,~ 9„9,„, -\i/2 



A(0) = 



(2TT 2 n 2 \U\0 D a 2 - 1) 



(16) 



(17) 



This self-consistent solution vanishes when a — ai — 
(2t: 2 Q 2 0d\U\)~ 1 ^ 2 , which is the minimal island area. Let 



us estimate the minimal island size, £ Q 



1/2 



for a re- 



alistic system, composed of a Nb island embedded on a 
semi-conducting quantum- well made of Si or GaAs. Tak- 
ing the effective mass m* jm = 0.98 and m* /m — 0.063 
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for Si and GaAs respectively* 2 ^ one can estimate the 2D- 
DOS in the quantum well. Taking for Nb T c = 9.26K 
and 6 D = one finds that for the Nb/Si hybrid 

the minimal island radius is £/ ~ 80nm, and for the 
Nb/GaAs system £j ~ 20nm. Both these lengths are still 
in the point-island regime, since the Fermi wave-length 
may be an order of magnitude larger for such quantum 
wells. 

Eq. (|17p also supplies us with a dependence of A on 
the island size and interaction strength. In the inset of 
Fig. [3] we plot a fit of the numerical data in the region 
U ~ U c (stars) to a function of the form A(?7) oc (U/U c — 
l) x O(U - U c ) (solid line), as in Eq. (17]). The fit yields 
the exponent x = 0.5002, in good agreement with the 
continuous model. 

In the definition of the model (Eq.©-©) we have ne- 
glected the boundary conditions on the normal-SC inter- 
face. Omitting the boundary effect is hard to justify a- 
priory, especially when the size of the SC grain is smaller 
than the SC coherence length. Accounting for the bound- 
ary should be accounted for by solving Eqs.([T5 ]) -(fT5 |) with 
an additional constraint A (a) = 0. However, the com- 
parison between the numerical calculation (in which the 
boundary condition are inherently implemented) and the 
analytical result (inset of Fig. [3]) shows a striking equiv- 
alence between them. This indicates that neglecting the 
boundary effect merely results in a quantitative modifica- 
tion. It does not change the qualitative behavior, which is 
mainly characterized by the power-law dependence spec- 
ified in Eq.([T7|). 



V. DISCUSSION 

One main feature of the result shown in Eq. (jTTJ) is 
that, in contrast to previous works on the proximity 
effectiiSiii, the length-scale is not the usual supercon- 
ducting coherence length £, but rather a new length scale 
ai = (2tt 2 £1 2 0d\U\)~ 1 / 2 . In order to understand the ori- 
gin of this new length scale we cast the criterion for 
the vanishing of SC correlations in the island given by 
Eq. JT7D to the form 



±(\E1 

V2\0D 



1/2 



2na 2 fl0D = f 



(18) 



The factor 2-kcl 2 Q,0d is nothing but the number of elec- 
trons with energy in the range d within the island, and 
thus the condition turns out to be a continuous version of 
the Anderson criteria, which states that SC vanishes in 
a grain once the number of pairs, roughly given by A/6 
where S is the level spacing, becomes less then unity. 
However, in the above model there is no discreteness of 
energy levels. Rather, the number of pairs in restricted 
due to the finite region to which the interaction is limited 
to. It is also clear from this argument why £ does not 
play a role in this system, as £ indicates the existence of 



a region where superconductivity is developed to its bulk 
value, which is not the case here. 

Yet another way of understanding the result of Eq. (fT7|) 
is to note that in a SC-semi-conductor junction, one ex- 
pects that due to the low density on the normal side, 
the suppression of pair correlations in the SC (due to the 
proximity effect) will no longer be on a length scale £ as 
in a SC-metal junction, but rather a new length scale, 
which in the point-island approximation corresponds to 

1 /2 

£r ps a j . Thus, if the SC island is smaller than the 
length-scale on which pair correlations are suppressed, 
SC in the island will vanish. One also expects that upon 
increasing the island size (beyond the point-island ap- 
proximation) this new length-scale will change. This 
problem is beyond the scope of the present work and 
will be addressed in a future study. 

In Eq. [TTJ A has a power-law dependence on the 
grain size. This is in contrast to the case of an iso- 
lated grain, where an exponential dependence on size is 
predicted^. Furthermore, due to the re-normalization 
of electron number and the dependence of critical island 
size on the DOS, the critical size may be either larger 
or smaller than given by Anderson's criteria. This may 
affect the possibility of fabricating devices made from SC 
grains embedded on a metallic matrix. This effect may 
be tested experimentally by, e.g., varying the DOS of 
the metallic substrate by changing its carrier density (by 
gating the sample, for instance). 

More intuition on the existence of a critical island size 
or interaction strength may be obtained by noting the 
similarity between Eq. [XT] and Anderson's criteria for the 
existence of a magnetic impurity in a metal^ 2 .. This sim- 
ilarity, along with the Friedel-like oscillations of Eq. [T2l 
implies that the Negative-U impurity is screened by the 
free electron gas, in an analogous way to the screening 
of a magnetic impurity. It would thus be intriguing to 
investigate the possibility of the formation of an effect 
equivalent to the charge Kondo effect—, resulting from 
the presence of embedded SC grains. 

We conclude by noting that these results may be tested 
experimentally, by planting superconducting impurities 
on a metallic substrate and measuring the local gap, in 
a similar way to that of Ref. U Another system in which 
our results may be valid is a SC grain strongly coupled 
to matching leads (i.e. SC Aluminium grain and nor- 
mal Aluminium leads), where the existence of SC may 
be verified as a function of grain size. While this may 
be experimentally challenging, it may be achieved by, for 
instance, planting magnetic impurities in the leads. The 
dependence of the SC gap on the properties of the sub- 
strate, as seen in Fig. [5] for instance, may be tested by 
changing the density on the metallic substrate. 
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Meir and Y. Avishai, and is grateful to T. Aono for help- 
ful discussions and for carefully reading the manuscript. 
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